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This note shows that Heisenberg's choice for a wave function in his original paper 
on the uncertainty principle is simply a renormalized characteristic function of a 
stable distribution with certain restrictions on the parameters. Relaxing Heisenberg's 
restrictions leads to a more general formulation of the uncertainty principle. This 
reformulation shows quantum uncertainty can exist at a macroscopic level. These 
modifications also give rise to a new form of Schrddinger's wave equation as the 
equation of a vibrating string. Although a heat equation version can also be given, 
the latter shows the traditional formulation of Schrddinger's equation involves a 
hidden Cauchy amplitude assumption. 
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A generalized wave packet 

We begin by showing that Heisenberg's choice for a wave function in his original paper 
^| on the uncertainty principle is simply a renormalized characteristic function of a stable 
distribution, S a ^(x; m, c) with a = 2 and (3 = 0, and location and scale parameters m and 
c. Relaxing the assumptions on a, (3 so that < a < 2,(3 ^ 0, leads to a more general 
formulation of the uncertainty principle. These modifications also give rise to a new form 
of Schrddinger's partial differential equation. 

Consider the following wave packet ip(x, t), where at time t = 0, ip(x, 0) has the form 

ip(x, 0) = A a exp[imx — c\x\ a ], (1) 

where 

It is easy to see that ip(x, 0) is normalized to unity: 

/oo roo roo 

* 0, 0)i/j(x, 0)dx = A 2 / exp[-2c\x\ a ]dx = 2A 2 a / exp[-2cx a ]dx. (3) 
-oo J — oo JO 
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Using the relation 

r y k e-y a dy = ±r(ttl) (4) 
jo a a 

and making the substitution u = (2c) *x, we obtain 

poo 2A 2 r°° 2A 2 1 1 

2A 2 a / exp[-2cx a ]dx = °- / e~ u du = ^~r(-) = 1. (5) 

°Jo FL J (2c)- Jo (2c)i« V 

Now, the form of the wave packet in Eq.(l) can be compared to Heisenberg's original wave 
packet, denoted here H(x, 0): 

H(x, 0) = (2r)3 exp[2nia x - ittx 2 }. (6) 

If we make the substitutions 

2na = m (7) 

7TT = C (8) 

a = 2 (9) 

in -^(x, 0), the wave packet of Eq.(l), we obtain H(x,0). (Note that with a = 2, A Q = 
[ff ]i = (f )i = (2r)l .) 

Now let's derive the amplitude function of ^(^ 0), which will necessarily also give us the 
amplitude function of H(x,0). First note that the log characteristic function of a stable 
distribution is 

/oo r£ iprt 
exp[ixz]dF( ) (10) 
-oo C 

= imz - |c'| Q |z| a [l + i/3(z/|z|)tan(7ra/2)], if a ^ 1 (11) 
= imz - \c\ a \z\ a [l + i/3{z/\z\){2/it)\og\z\], if a = 1 (12) 

where m is a real number, d > 0, < a < 2, 1/31 < 1. Proof of this theorem, due to 
K^^ne and Lev y in 19 36, ^ be found in"Q „ r Q _ He re a, the cU^Uc 
is essentially an index of peakedneess (a = 2 for the normal or Gaussian distribution, 
a = 1 for the Cauchy distribution). The parameter (5 is an index of skewedness {(3 = 
for symmetric distributions). The parameter d = c» is a scale parameter (the standard 
deviation when a = 2). Finally, m is a location parameter (the mean if a > 1; it is also the 
median or modal value of the distribution if /3 — 0). 
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For (3 = we obtain the characteristic function of a symmetric stable distribution, which 
is identical to Eq.(l), if we omit the normalizing constant A Q . Therefore, for the amplitude 
function of our wave packet, we take the Fourier transform, A(z), of Eq.(l) to obtain 

/oo rco 
tp(x, 0) exp[— ixzjdx — A a exp[imx — c\x\ a ] exp[— ixz]dx = A s a fi(z; m, c). 
-CO J — CO 

(13) 

In other words, we obtain a symmetric stable density function s a fi(z; m, c) = dS a fi{z;m,c) 
with the normalization constant A for the amplitude function. The symmetric stable density 
has < a < 2, (3 = 0, and location and scale paramters m and c, respectively. Note that 
the amplitude function is normalized so that the integral of its square is equal to 1. This 
involves the square of the probability density function s a fi(z; m, c). 

For Heisenberg's case where a = 2, we may explicity solve for A(z) = A(a), which will 
be necessarily the Gaussian density multiplied by a normalizing constant. We reintroduce 
a factor of 2n to obtain 

/co 1 2 i 7l(<J — CT ) 2 
(2r)4 exp[27ricr a; — nrx 2 ] exp[— 2nixa]dx = (— )* exp[ — ]. (14) 
-oo T T 

This is Heisenberg's amplitude function. That the integral of its square is 1 follows from: 
/ [A(a)] 2 dar= (-)2 / exp[ K — ; -}da = (-) 2 (^- ) 2 2 / e du = ^= = 1, 

J — co T J — co T T Z7T JO V 

(15) 



where we have used the substitution u = y ^r(cr — c ) . Note that the usual normalizing 
constant for the Gaussian distribution (where d is the standard deviation) has been 
absorbed into A . So above and below, when we write the stable density s a ^{z\ m, c), 
we will understand the omission of the usual normalizing constant, and will consider only 
the normalizing A Q in the product A s a ^(z;m,c). This will ensure that the square of the 
amplitude function is a probability distribution. 

Alternative amplitude functions 

For a — 1, which corresponds to the Cauchy distribution, the normalizing constant 
A Q = c5 = ((/)? = (c / )a, so the amplitude function is 

i I~2 c 

A(z) = A o s 1:0 (z;m,c) = c*\l — ^— rx, (16) 

V 7r & + [z — my 
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where we have removed a division by \j2n in the usual statement of the Cauchy. That this 
is the correct normalization for the amplitude function in Eq.(16) follows from the integral: 

2 f°° 2 f°° I 

} A{zfdZ = n L [c^{z-mYY dZ = n L jjTir*"' ^ 

where we have used the substitution y = ^p 1 - We may now appeal to the relations, for 
a, c > and n a positive integer: 

r dx 1 x 2n — 3 r dx 

J (ax 2 + c) n ~ 2(n - Tjc (ax 2 + c)"- 1 + 2(n - l)c J (ax 2 + c)^ 1 ^ ' 

and 

/_^_= 1 tan- 1 ^]. (19) 
J ax 2 + c vac V c 

Thus we get 

2 r°° dy 2 1 r°° dy 2 1 r°° dy 2 , _ . . 

--2/ — ^ = - tan" 1 ^ = 1. (20) 



7T 7-oo [1 + y 2 ] 2 7T 2 7-oo 1 + y 2 n 2 Jo 1 + y 2 it 
If we further generalize Eq.(l) by relaxing the constraint on f3, we obtain the wave function 

ip(x,0) = A Q exp{[imx - c\x\ a ][l + i(3(z/\z\)tan(na/2)]}. (21) 

Note for the wave function in Eq.(21) that since % multiplies /3, the normalizing constant 
A Q given in Eq.(2) is unchanged in terms of a. For a — |, which we will now consider, 
A Q = c = (c')5. Thus for a — \ and (3 = — 1, we obtain for the amplitude function 
the completely positive stable distribution (sometimes called Pearson V), multiplied by the 
normalizing constant c: 

A(z) = As^z^c) = o ^ z C _ m)3 «p["2(7r^)]- ( 22 ) 

As a check, we integrate the probability function P(z) = [A(z)] 2 corresponding to the 
amplitude function in Eq.(22): 

= 2j( M V« 2 = 2-r(-) = 1, (24) 

where we have used the substitution u — r. 

(z-m)2 
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Finally, for the general case, we may express the amplitude function as a renormalized 
stable density, which is in turn represented by a Taylor expansion in the form of gamma 

n n 

functions j2j[p- 583] (alternative expansions may be found in 



A(z) = A oSat(3 (z; 0, 1) = A -y- E fc , (-*)* sm [^(^ - »)]> 



(25) 



for z > and 1 < a < 2. For z < we have the general relation s a>/ g(— 2; m, c) = 
s a -^(z; m, c). For < a < 1 we have the similar expansion, for z > 0, 

A(*) = A oSqi/3 (z; 0, 1) = A -^f / T{1 ^ ka \ -z- a ) k sinfo/3 - a)]. (26) 

z V 7T fc=1 ft! 2 

We may recover m and c in Eqs. (25,26) by the substitution z — - \ l . 



The uncertainty relation 

Now let's consider the uncertainty relation. From Eq.(l), where the distribution is sym- 
metric, we get the value for (Ax) 2 as: 



(Ax) 2 = / i)*(x,0)x 2 ip(x,0)dx. (27) 
Inserting a factor of u 2 = (2c) «x 2 into the calculation of Eq.(5), we obtain 

1 r(^) 

(Ax) 2 = 2 ^ ? . (28) 

For a = 2 this yields (Ax) 2 = 7-, or in Heisenberg's formulation . 

Next consider the uncertainty in z (or a). First consider the case a = 2. From Eq.(15) 
we have 

(Aa) 2 = T (a - a ) 2 [A(a)} 2 da = f° (a - <x D ) 2 exp[- 27r(t7 " = f • (29) 

J-00 T J-00 T 47T 

Thus we obtain the uncertainty relation 

AxAa = (30) 

47T 

From the de Broglie relation Ap = hAa, where h is Planck's constant, this becomes 

AxAp=^. (31) 
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However, for comparison with the results below, we will use for the (renormalized) Gaussian 
amplitude, the uncertainty relation in the form 

AxAz = 1. (32) 

Note that for the Cauchy density, where a — 1, f3 — 0, the mean and variance don't 
exist ("are infinite"). But we are considering a Cauchy amplitude, and hence the square of 
the Cauchy density (renormalized) for the probability density. For this density the second 
moment exists, as we will now demonstrate. From Eqs.(16,17), we calculate (Az) 2 as: 

(A,)* = £> - nOW* = I £ [c ,%_%^ (33) 

2c 2 Z" 00 v 2 , 2c 2 1 r°° dy 9 

dy = —-2 / — y — = c 2 , (34) 



7T J-oo [1 + y 2 ] 2 7T 2 Jo 1 + y 
where we have used the relation 

x 2 dx 1 x 1 f dx 



(35) 



(arc 2 + c) n 2(n - l)a (ax 2 + c)™" 1 2(n - l)a 7 (ax 2 + c)"- 1 
Thus we obtain the uncertainty relation, from Eqs. (28,34), 

AxAz = -^=. (36) 

For the Pearson V amplitude, we have from Eqs. (22,23) 

{z-m) 2 [A(z)] 2 dz = c 2 j— — - e - — dz = c 4 -e~ y dy, (37) 
-oo Jm yz — in) Jo y 

where we have used the substitution y = ^ z c _ m ^ ■ This integral is divergent. So instead we 
calculate 

/oo roo f£ c 2 /■oo 

\z-m\[A(z)] 2 dz = J — — —e~?^dz = c 2 J^ e' y dy = c 2 . (38) 

This yields, from Eqs. (28,38) the uncertainty relation 

AxAz = (39) 
It is easy to see from Eq.(28) that the general uncertainty relation, as a function of a, is 

AxAz 



i rm 

{a) (40) 



N(2)= r (=)" 

This, then, is the reformulation of Heisenberg's uncertainty relation. The uncertainty is a 
function of the characteristic exponent a of the (renormalized) stable amplitude. As a — > 0, 
the uncertainty becomes unbounded. 
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The time-dependent wave function and the dispersion relation 

We can write the time-dependent wave equation corresponding to Eq.(l) as a superposi- 
tion of plane waves: 

/oo 
A(z) exp[i(zx - v{z)t)]dz, (41) 
-oo 

where A(z) is the stable amplitude — a renormalized stable density, and v{z) is the frequency. 
A dispersion relation connects v(z) to z. 
From the de Broglie relations 

E = hu (42) 
p = hz (43) 

we obtain the relation 

u = z—, 44 
V 

which gives as the time-dependent wave equation 

f°° E 
ip(x,t) — / A(z) exp[iz(x t)]dz. (45) 

J-oo p 

2 

(Note that we do not insert the classical relation E = j^, where M is mass, at this point, 
because doing so does not yield a proper inverse Fourier transform.) Each plane wave 
equation g(x,t) = exp[iz(x — ^t)\ has differential operators ^ and with eigenvalues 
— z 2 and — z 2 ^- respectively: 

d 2 9 



dx 2 

&g = _ z2 E 2 
dt 2 p 



-z 9 (46) 



= -z 2 -g. (47) 



These relations give rise to the partial differential equation 

d 2 g E 2 d 2 g 



(48) 



dt 2 p 2 dx 2 

The time-dependent wave equation in Eq.(45) may be rewritten more fully (for a ^ 1) 



as 



EE EE 

Mx,t) = A exp{[im{x t)-c\(x t)\ a ][l + i(3((x t)/\(x t)|)tan(7ra/2)]}. (49) 

p p p p 
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For symmetric distributions (f3 = 0), the probability density function corresponding to 
ip(x, t) is 

P{x,t) = ip*{x,t)ip(x,t) = A 2 Q exp[-2c|(ar- jt)\ a ], (50) 

which is the characteristic function of a stable density. For a — 2, this is the Heisenberg 
time-dependent density. 

Schrodinger's equation revisted 

Schrodinger's equation may be viewed as a simple consequence of the Heisenberg uncer- 
tainty relations. Eq.(49) is a solution of the partial differential equation Eq.(48), so we have, 
as replacement for the Schrodinger equation, the partial differential equation 

d 2 ip E 2 d 2 ip 



dt 2 p 2 dx 2 ' 

which may be rewritten in the form 



(51) 



This, of course, is the equation of a vibrating string, where v = - is the speed of propa- 



v 



gation of the waves. It is a true wave equation, by contrast to Schrodinger's heat equation 

dt dx 2 



formalism, which relates ^ to In fact, noting from Eq.(49), letting (3 equal zero for 



simplicity, and letting sgn y denote sgn (x — —t) , that 



— — = (im — calx sgn y)ip (53) 

ox V 

= ff im _ ca \ x _ —t\ a ~ l sgn yf - ca(a - l)\x - -t\ a ~ 2 )i) (54) 

ox P V 

dip E dip 

dt p dx 

d 2 ip E 2 d 2 ip 



(55) 
(56) 



dt 2 p 2 dx 2 

it does not appear to be particularly useful to relate ^ to although this can be done 
In fact, 
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Only in the case of the Cauchy amplitude a = 1 do we find this latter formulation in a 
simplified form: 

d± = E 1 ^ (58) 
dt p (im — c) dx 2 

2 

If we now make the substitutions E = ^g, p = ha we obtain 

dip ha im + c d 2 ip 



dt 2M (m 2 + c 2 ) dx 2 

which may be rewritten 



(59) 



ch/> h 2 a m — ic d 2 ip 
~dt = ~2M{m 2 + c 2 ) dx 2 ' ( ' 

It would appear that the traditional Schrodinger equation involves a hidden Cauchy ampli- 
tude assumption. The latter equation can be divided into two equations, one involving m 
and the other involving — ic. 



Conclusion 



Stable distributions are the only distributions that exist as limit distributions of sums of 

random variables, thus giving rise to central limit theorems. Therfore they play a paramount 

role in the physical world. We have shown that Heisenberg's original choice for a wave 

packet to illustrate his uncertainty principle is simply the characteristic function (the inverse 

Fourier transform) of a Gaussian distribution, leading to a Gaussian amplitude function 

with a = 2 and f3 — 0. Relaxing Heisenberg's assumptions to the general case < a < 

2, \(3\ < 1, leads to stable amplitudes renormalized so that the integral of their squares 

are probability distributions. The renormalization constant gives rise to a new form of 

Heisenberg's uncertainty relation, expressed in terms of the characteristic exponent a of 

/ 1 r(— ) 

the underlying stable amplitude: AxAz = J — TT7T\- This relationship was illustrated 

y (2) a 1 ( a ) 

by explict calculation for the Gaussian (a = 2), the Cauchy (a = 1), and the Pearson V 
(a — |, j3 — —1). As a — > 0, the uncertainty AxAz becomes unbounded. This means 
that, depending on the underlying stable amplitude, quantum uncertainty can arise at a 
macroscopic level. 

2 

By eschewing the ad hoc classical insertion E = ^g, we were able to solve for the time- 
dependent wave equation as a superposition of plane waves, by taking the inverse Fourier 
transform of the stable amplitude function. For a = 2, this recovers Heisenberg's case. The 
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wave function follows the partial differential equation = ^2 |^r, which is the equation for 
a vibrating string. This is a proper wave equation, differing from Schrodinger's equation, 
which is really a heat equation as it relates instead of Jp-, to -£-5. The traditional form 
of the Schrddinger equation can be recovered, but only in the case a = 1. Thus it would 
appear that Schrodinger's equation involves a hidden Caucy amplitude assumption. This is 
not fatal, but is limiting. The more general heat equation relationship is given by Eq.(57). 



Email: quantum@orlingrabbe.com 



[1] Bergstrom, H., 'On some expansions of stable distribution functions', Arkiv fur Mathematik, 2, 
1952, 375-378. 

[2] Feller William, An introduction to probability theory and its applications, Volume II, 2nd edi- 
tion, John Wiley & Sons: New York, 1971. 

[3] Gnedenko B.V., and A.N. Kolmogorov, Limit distributions for sums of independent random 
variables, Addison Wesley, 1954. 

[4] Heisenberg W., 'Uber den anschaulichen Inhalt der quantentheoretischen Kinematik und 
Mechanik,' Zeitschrift fur Physik 43, 172-198 (1927); reprinted in Dokumente der Naturwis- 
senshaft, Vol 4 (1963), pp. 9-35. 

[5] Levy, Paul, Theorie de V addition des variables aleatoirs, Gauthier-Villars, 1937. 



